Second magnetization peak in flux lattices - the decoupling scenario 
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The second peak phenomena of flux lattices in layered superconductors is described in terms of 
a disorder induced layer decoupling transition. For weak disorder the tilt mudulus undergoes an 
apparent discontinuity which leads to an enhanced critical current and reduced domain size in the 
decoupled phase. The Josephson plasma frequency is reduced by decoupling and by Josephson glass 
pinning; in the liquid phase it varies as 1/[BT(T + To)] where T is temperature, B is field and To 
is the disorder dependent temperature of the multicritical point. 



Vortex matter in the presence of disorder has emerged as a fundamental problem of elastic manifolds in a random 
media (!]]. Impurity disorder does not allow long range translational order of the flux lattice and finite domains are 
expected At low temperatures and fields the system is a Bragg glass i.e. the lattice is dislocation free, 

at long scales the displacement correlations decay as a power law and Bragg peaks are expected. The impurity 
induced domains are essential for the description of both equilibrium, e.g. thermodynamic phase transitions and 
non-equilibrium, e.g. critical current phenomena. 

The critical current j c measures the pinning force in the domains Jl],|| . Increasing the magnetic field or temperature 
reduces the pinning force and j c is decreased. However, in many type II superconductors a sharp enhancement of j c is 
observed at a "second peak" field Bo- This peak phenomena is most pronounced in layered superconductors such as 
Bi 2 Sr 2 CaCu 2 0s (BSSCO) §-§, YBa 2 Cu 3 7 (YBCO) @, in NbSe 2 @[lC]] and in Pb/Ge multilayers (ll) for fields 
perpendicular to the layers. The second peak phenomena signals that pinning becomes more effective, e.g. due to 
softening of the flux lattice 01. The reason for softening could be the approach to melting fl2| ; however, neutron 
scattering data on BSCCO j^j shows that Bragg spots of the flux lattice persist well above Bq. 

Disorder plays an essential role also in the equilibrium phase diagram of layered superconductors. This has been 
most extensively studied in BSSCO l5|~[^ JT^ , p^|| . The second peak corresponds to a phase transition [Q in the range 
500 — 900G (decreasing with disorder) and is weakly temperature dependent up to a temperature To iOK. The 
point Bq, To is a multicritical point where the second peak transition meets a first order transition as well as two 
depinning lines. Thus the second peak manifests both equilibrium and nonequilibrium phenomena of disorder in flux 
lattices and its understanding presents a fundamental challenge. 

For a flux lattice with point impurities, by using renormalization group (RG) and replica symmetry breaking (RSB) 
methods we have derived |l6| a phase diagram with four phases, which all meet at a multicritical point Bo, To, in 
remarkable correspondence with data on BSCCO. The present work focuses on the layer decoupling transition at a 
temperature independent field Bo for T < Tq. As shown here, the fusion of Bragg glass concepts with decoupling 
accounts for the peculiar second peak phenomena, i.e. the enhanced j c . The Josephson plasma resonance is also 
considered as a probe of the Josephson coupling |L7|,[l8), being reduced by decoupling and by a Josephson glass 
parameter. Very recent data on BSCCO has indeed shown a significant reduction in the resonance frequency at the 
second peak transition . 

It has been recently shown that decoupling coalesces with a defect unbinding transition pl| which has analogs in 
isotropic systems j2^]. The resulting vacancies and interstitials lead to a reduction in the elastic tilt mudulus p3[ , 
consistent with the decoupling scenario as described below. It is possible then that a decoupling-defect transition 
accounts for the peak phenomena in all type II superconductors. The analysis below is, however, presented for layered 
anisotropic systems where quantitative predictions can be made. 

In a layered superconductor each flux line is composed of one point singularity, or a pancake vortex, in each layer. 
When the pancake vortices fluctuate they can generate a divergence in the Josephson phase, leading to a renormalized 



Josephson coupling Ef which vanishes in the decoupled B > Bo phase |16|,^4 25 1. The 3-dimensional flux lattice is 
still present in the decoupled phase (in the Bragg glass sense), with interlayer coupling mediated by the magnetic 
field. Before presenting a microscopic model, I start with a rather simple description of elasticity within domains, 
which shows the second peak transition, i.e. j c enhancement at decoupling. 

The transverse tilt modulus of a flux lattice in a layered superconductor for fields perpendicular to the layers is 
given by |f|§ 

C44(q < k) = 32^d + S l + W + A^ + (S^ ln(fl2/4< ° 2) (1) 
where q and k are momenta parallel and perpendicular to the layers, respectively, X a b and A c are the London 
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penetration lengths parallel and perpendicular to the layers, respectively, (j>o is the flux quantum, a 2 — <fio/B is the 
unit cell area, d is the interlayer spacing, £ is the in-layer coherence length and r = (j^d / {Air 2 }? ab ) sets the energy 
scale. The first term of Eq. (|f|) is due to the magnetic coupling, while the 2nd and 3rd terms originate from the 
Josephson coupling energy per unit area Ej, i.e. X 2 — r\ 2 b / (AnEjd 2 ). The second term is peculiar: at q ^ it 
vanishes when Ej vanishes and A c oo, as it should. However, at q = this term seems to survive even if A c — > oo. 
The origin of this peculiarity is that the harmonic expansion of the Josephson cosine term which identifies C44 fails |2q] 
when both q, 1/A C — > 0. In fact, the nonlinear cosine term generates a renormalized A,? which diverges at decoupling. 

The Bragg glass domain size Rbg (parallel to the layers) sets a scale for the relevant q values. When Rbg > A^ 
the tilt mudulus is large, containing the B 2 /Air term of Eq. (|l|). However, as decoupling at the field Bq is approached 
A,? diverges and when Rbg < A,? Eq. ([!]) fails to describe C44 on the scale of q sa \/R B g- This defines an anharmonic 
crossover regime where usual elasticity cannot be used to derive Bragg glass properties. Finally, at B > B elasticity 
is restored and C44 is reduced to the first term in Eq. (|l|). The main interest is in the regime of strong fields, i.e. 
a < 2A a fc where To < t is below melting |16[ ]. Thus at B < Bo and for sufficiently large domains the second term in 
Eq. (1) dominates and C44 has an apparent discontinuity, 

c 44 = TrA^r/da 4 Af < Rbg (2a) 

c 44 = T/(32ir\ 2 ab d) Af = 00 (2b) 

Hence C44 is reduced within the anharmonic regime by the small factor e = a 4 /(327r 2 A 4 b ). 

The apparent discontinuity in C44 affects also the domain sizes which can be estimated by a dimensional argument 
[§,U ■ Consider the tilt C44 and shear cqq terms of the elasticity Hamiltonian for the displacement u(r) and its transverse 
component u^(r). Rescaling parallel and perpendicular lengths yields an isotropic form plH 



H = f dMKf c 66 3 [VMr)] 2 - (eo/a 2 d)U pm (r)J2™sQ ■ [p-u(r)]} (3) 

J Q 

where U P i n (r) is a random potential in 3-dimensional r = (p, z) which couples to the flux density modulations with 
wavevectors Q; its disorder average is (U P i n (r)U P i n (r')) — ^dUS 3 (r — r'). Disorder average over configurations u(r) 
and u'(r) yields ~^2q cosQ • [u(r) — u'(r)]; the sum is cutoff by Q < (m 2 ^) -1 / 2 where (u 2 ) w (u^) are the fluctuations 
in a domain of size R'. Thus averaging Eq. (0) yields 

(H)/R> 3 = \ c ^ c %\u 2 T )R'- 2 U^ 2 e /[a 2 d(u 2 T )R^ /2 ■ (4) 

Minimizing with respect to R' yields R' ~ (w 2 -) 3 , i.e. the Flory exponent ||. The domain size parallel to the layers 
is (up to hi(a/d) and a numerical prefactor) 

R^(X ab /a) 5 (u 2 T ) 3 /(s^d) \?<R 

R « (A afc / a ) 3 (4) 3 /(47r^ 4 d) Af = co (5) 

where c 66 = r/Wda 2 p6|-p8t, s = AnpU\\ b /[T 2 a 2 ln 2 (a/d)] defines the decoupling transition at s = | and p w I is 
defined below. 

The pinning length R = R p is given by Eq. (|) with {u 2 T ) « £ 2 . To allow for large pinning domains one needs either 
a < A a b or to allow for domains with a somewhat larger fluctuations in (u 2 -,); the latter increases R p very rapidly since 
it increases with the 6-th power of ut- The critical current can now be estimated by balancing the Lorenz force 
j c BR 3 jc with the pinning force (H)/£o (evaluated at the minimum of Eq. (^)), leading to j c ~ I/C44. Increasing the 
field within the anharmonic regime decreases C44 by the factor e so that j c is significantly enhanced when a < A a ft. 
Note that the measured magnetization changes (and inferred j c ) at Bq decrease with temperature due to the strongly 
temperature dependent relaxation rates approaching the much smaller equilibrium magnetizations. 

A second length scale R = Rbg is identified by Eq. (||) with the fluctuations (u^) « a 2 . The proper definition of 
Rbg is the scale for the onset of the lnr form for the displacement correlation function. While the derivation from Eq. 
(||) cannot capture this lnr, it does give the right expression for Rbg J|]- Thus, Rbg depends on C44 and is reduced 
by e 1 / 2 through the anharmonic regime. The latter depends also on A^for which In A,? ~ (B — Bq)^ 1 in the RSB or 
1st order RG solutions [jlfj), though In A,? ~ (B — Bq)^ 1 / 2 in second order RG [Q; decoupling may also be of first 
order p5{| , leading to a narrower anharmonic regime. Fig. 1 illustrates the lengths Rbg an d A^, demonstrating the 
anharmonic regime within which Rbg A a s a significant drop and correspondingly j c has an apparent jump. Note that 
even in the decoupled phase (B > Bq) Rbg is large for typical type II superconductors, Rbg ~ ^ab a3 /(4 7rs £o c ^ a ' 
consistent with a decoupling transition within the Bragg glass phase, i.e. below a melting transition. 
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FIG. 1. Bragg glass domain size Rbg parallel to the layers and the renormalized London length perpendicular to the layers 
A,?; the latter diverges at the decoupling field Bq. Rbg can be found from elasticity for B < Bo only if Rbg > A,?; otherwise, 
as in the hatched region, the elastic tilt mudulus is ill defined. 

I proceed now to derive the lattice displacement correlation allowing for a renormalized Josephson coupling and 
for a Josephson glass order parameter. This derivation avoids the harmonic expansion for the elastic modulii and 
shows how the Bragg glass domain sizes are directly affected by the renormalized A,?. The Josephson phase between 
the layers n and n + 1 at position r in the layer involves contributions from a nonsingular component 9 n (r) and from 
singular vortex terms ]30| ] . Consider a flux lattice with an equilibrium position of the Z-th flux line at R; . The singular 
phase around a pancake vortex at position R/ + u™ is a(r — R/ — u") where a(r) = arctan(y/ai) with r = (x,y). 
Expansion of the interlayer phase difference a(r — R; — u") — a(r — R/ u ™ +1 ) yields for the singular part of the 
Josephson phase b n (r) = X^( u ™ +1 — uf)Va(r — Hi). The Hamiltonian for the transverse displacements involves also 
the magnetic contributions to the shear modulus cqq = r/(16da 2 ) and the tilt mudulus 



Here c(q, k) = (a 2 /27rd) 2 [k 2 c1 4: (k) + q 2 cee]/Tk 2 , Ej is the Josephson coupling and the coefficient of the non-singular 
phase is @ G f (q,k) = <±ird?T(\- 2 + k 2 )/(Tq 2 ). The conventional C44 is obtained by expanding the cosine term in 
Eq. (|^) and shifting 9(q 1 k) to eliminate the cross term. The latter shift leads to an expansion parameter |2^| with 
terms ~ q 2 k 2 \uT(ci, k)\ 2 /[q 2 + A,T 2 (1 + A^fc 2 )] 2 , i.e. these diverge when both q, 1/A C — > and the expansion becomes 
invalid. 

Consider now a pinning potential Up in (r) which couples to the vortex shape function p(r) leading to a pinning 
energy J d 2 r J2 n 1 Up in {v)p(r — R; — u") . The aim is to identify domain sizes R p (and infer Rbg)i hence the 
pinning energy is expanded in u™ and a replica average with the weight exp{— J d 2 rJ2 n [Upi n ( r )] 2 /U} then leads to 
exp[(f7p/4T 2 ) l g u" ,a • u™'^] where J dip(r)djp(r)d 2 r = pS^j and a,fl= 1, 2, . . . , n are replica indices. 



The b a (q,k) variables can be decoupled from the total Josephson phase b n (r) = b n (r) + 9 n (r) by shifting to 
d a (q,k) = b a (c[,k) - B 7!a (q,/c)G7 1 (q,fc)6 tt (q,fc) where 




S aj?( q ' k "> = G f 1 ^ k ) a (^ k)5 a ^ - s q 2 /k 2 z 



a(q, k) = 1 + G/(q, fe)c(q, k)q 



' 2 and so = Upa 2 d/ {And 2 T) 2 . The resulting replicated Hamiltonian is 
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H r = \ ]T B-^d a (q,k)dP%q,k) + ±{ C (q,k)a-\q,k)qH a , - S oa-\q,k)q 2 ^ 



q,k;ct,(3 

" ^ E / d2 rcosK(r) - ^ £ / d 2 r cosfe(r) - S£(r)] (7) 

The inter-replica i?„ term is generated from the Josephson coupling in second order RG. It is essential to keep it from 
the start since it generates a Josephson glass parameter and affects the value of the decoupling field Jlj| . 

The a(q, k) factor, which results from the nonsingular phase, is for eCl very close to 1 for all q, k values except 
when both k < 1/X a b and q > fca/A a {,. The phase transitions are dominated by k > l/a modes so that our previous 
phase diagram is recovered ( pM Fig.l). In particular there is a multicritical point at a field Bq where s = | 
and temperature T = to 2 ln(a/c?)/87rA 2 b . At B = B and T < T there is a decoupling transition at which the 
renormalized Josephson coupling z (with bare value Zb are = Ej/Td) vanishes. Note that the higher Bo of YBCO as 
compared to BSCCO is consistent with a shorter A t and a somewhat weaker disorder. 

The fluctuations in Uy(q, k) in terms of the shifted variables, using the RSB solution jl6| are given by 

<Mq,fc)| 2 ) = (2nd 2 )- 2 ii^q 2 G f (c l ,k)a- 1 (c l ,k) ( c(q, fc) 9 2 + ^ ) 

k% y Gj {q,k) + zj 

q 2 + z) } + ... (8) 



s ( c(q, k) J, 



c(q, k)a 2 (q, k) \ a(q, k) 

where ... stands for terms which converge in (q, k) integration. Note the term Gj 1 (q, k)z/[Gj 1 (q, k) + z] which 
depends on the order of q — > and z — > limits; this limit dependence leads to the apparent discontinuity in C44 as 
discussed above. For z ^ and small q, i.e. Gj 1 (q, k) -C z the first term in Eq. (|J) dominates, leading to 

^ fafc)|a)W a-[c^W]' 9<1/Af (9) 

where C44 is from Eq. ( pa] ) and the condition Gj 1 (q, fe) <C 2 is written in terms of a renormalized London length 
\f = [\ 2 b r/(4:irTd 3 z)] 1 / 2 . The correlations at distance r parallel to the layers are then 

(Mr) - u T (0)} 2 ) « - ■ (10) 

a c 44 c 66 ^ 

The last equality defines the pinning length R p where the fluctuations become of order £p. This result for R p (up to 
a numerical prefactor) is the same as the one obtained from Eq. (||) with (v,j,) ~ a 2 . The Bragg glass domain size is 
enhanced by Rbg ~ ^p(a/Co) 6 , as discussed above. 

In the decoupled phase with z = the second term in Eq. (M) dominates. To leading order in e the result is identical 
to Eq. ( [To|) except that C44 is replaced by its z = value Eq. (ph|), i.e. the pinning and Bragg glass lengths are reduced. 
The main result is then that the fluctuations in ut{t) behave with an effective C44 which is large when q < 1/A^ (Eq. 
(pa|)), i.e. for domain sizes Rbg > A^, while for z — C44 is reduced (Eq. (pb|)). In the anharmonic region below 
decoupling (see Fig. 1) where Rbg < the full form of Eq. <JsJ> is required to interpolate between these limits; this 
form avoids the ill-defined harmonic expansion in this regime. 

Consider next the Josephson plasma frequency, given by lo 2 ^ — (c 2 /eoX 2 ) (cos b n (r) ) where eo is the dielectric constant 

|L7],[l8|]. The average in (cos 6„(r)) is on both thermal fluctuations and disorder and can yield significant information 
on the phase diagram. As shown by Koshelev jL8] the local (cos6„(r)) is finite even at high temperatures, e.g. 
above the decoupling transition. A high temperature expansion yields jl8| (cos6„(r)) = (Ej/2T) J d 2 r exp[— A(r)] 
where A(r) — Y^qkO- ~ cosc l ' r )(l^ a (q: k)\ 2 ). The solution with disorder |l6| yields (up to a ln_B dependence) 

A(r) = B(T + To)q 2 r 2 / {2BqTq) for r < l/q u where q u = 2ln 1/2 {a/d)/X ab while A(r) - lng„r or - r for larger r. The 
r integration is dominated by the short r correlation which yields 

{C ° Sbn{r)) ^2l^Jd)- BT(T + T ) 

A 1/ BT dependence has been obtained by Koshelev with a weakly temperature dependent prefactor for an XY 
model, i.e. infinite X a b- Data on BSCCO O] has shown that (cos6 n (r)) ~ jj—o.Srp—i j n reasona ble agreement with 
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the 1/BT form. The present result shows that in fact the l/BT form is valid in the disorder dominated regime, i.e 
T < T , though in general the fluctuation term yields u? vl ~ 1/[BT(T + T )]. 

Using the RSB solution, it can be shown that the Josephson glass parameter contributes a negative term to 
(cos6„(r)} so that lo p i is reduced, while the Josephson coupling contributes a positive w z/zb are term which vanishes 
at decoupling. These are mean field results to which fluctuation terms, as Eq. (|Tl|), should be added. The recent 
data on BSCCO |pj| , po| is consistent with these results, i.e. a drop at the second peak transition followed by a field 
dependent fluctuation term at higher fields. 

In conclusion, it is shown that a decoupling transition leads to an apparent reduction in C44 within an anharmonic 
region where the harmonic expansion fails. The proper interpolation across the anharmonic region is achieved by 
Eq. (||). The reduction in C44, the resulting reduction in domain sizes and the enhanced j c account for the hallmark 
feature of the second peak transition. Furthermore, B being weakly T dependent and decreasing with disorder ||, 
as well as the Josephson plasma resonance data [ p^|pOt |, lend substantial support for the identification of the second 
peak transition as a disorder induced decoupling. 
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